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Abstract 

The relation between the dilatation operator of M = 4 Yang-Mills theory and inte- 
grable spin chains makes it possible to compute the one-loop anomalous dimensions of all 
operators in the theory. In this paper we show how to apply the technology of integrable 
spin chains to the calculation of Yang-Mills correlation functions by expressing them in 
terms of matrix elements of spin operators on the corresponding spin chain. We illustrate 
this method with several examples in the SU(2) sector described by the XXXi/ 2 chain. 
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1. Introduction 

Recently there has been considerable progress in studying the AdS/CFT correspon- 
dence beyond the BPS and near-BPS limits (see [|J for a review). Integrability on both 
sides of the correspondence plays an important role in these quantitative checks. On the 
gauge theory side it has been shown that the planar one-loop dilatation operator can be 
identified with the Hamiltonian of an integrable spin chain ||. The vector space on each 
site of the spin chain depends on the particular sector one is interested in. For instance 
for scalar operators the dilatation operator can be described by the SO(6) spin chain in 
the vector representation 0], while an SU(2, 2|4) spin chain in the singleton representa- 
tion describes the full M = 4 Yang-Mills planar one-loop dilatation operator ||. The 
Bethe ansatz allows for finding (though implicitly) the exact planar one-loop anomalous 
dimensions of arbitrary local operators in gauge theory. 

On the string theory side of the correspondence, a different integrable structure has 
been found in the worldsheet theory [|J and was related to the gauge theory spin chain 
description in || . Nevertheless, finding the spectrum of string theory in AdS still remains a 
challenge. String theory becomes tractable in the semiclassical regime, when the quantum 
numbers of the string states become large The energies of certain semiclassical 

soliton solutions with two large angular momenta were shown to match the anomalous 
dimensions of operators with large quantum numbers computed to one loop using the 
XXX spin chain in || and to two loops using the elliptic spin chain in [Q. In this context, 
the higher local conserved charges on both sides have been matched in [HJ- Moreover, 



direct agreement has been shown between the continuum limit of the dilatation operator 
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in the coherent state basis and the Hamiltonian of the sigma model [1 l.| without needing 
the details of the semiclassical solutions. 

Less symmetric spin chains were also analyzed in relation to gauge theories. The XXZ 
chain with periodic boundary conditions was related to pure M = 2 SYM in |[12f. while the 



same chain but with twisted boundary conditions was related in [|13] to the (/-deformation 
of the M = 4 SYM. The same gauge theory was shown in [13 1 to be related to a limit of 
a multi-parametric spin chain which also describes certain non-supersymmetric theories. 
For related work on integrability in gauge and string theories see [fl5| - |5T |. 

So far these studies have been limited to finding the eigenvalues of the dilatation 
operator. In this paper we show how to use the technology of integrable spin chains to 
find three-point functions of arbitrary scalar operators in M = 4 Yang-Mills and some of 
its deformations. Even though it is relatively easy to compute three-point functions of 
operators at tree level, it becomes increasingly laborious to do so at higher loops. In fact, 
even the computation of the tree level three-point function coefficient of operators with 
definite one- loop scaling dimension is quite hard, because these operators are combinations, 
with occasionally rather complicated coefficients, of the natural tree level operators. The 
spin chain conveniently identifies a basis of operators with definite scaling dimension and 
leads to an algorithm for the computation of the leading three-point function coefficient. 
This coefficient turns out to be related to the expectation value of certain spin chain 
operators, a subject which has been extensively studied (see []59j and references therein). 
Loop corrections to the three-point function coefficient may be expressed as the products 
of expectation values of operators similar to the spin chain Hamiltonian. 

In §2 we construct the map between gauge theory and spin chain correlation functions, 
emphasizing the spin chain realization of arbitrary scalar operators. In §3 and §4 we review 
the relevant details of the algebraic Bethe ansatz and spin chain correlation functions. We 
proceed in §5 to describe several examples. 

At one loop the spin chain description of various sectors of the gauge theory is not 
restricted to operators with large quantum numbers. The planar perturbative correlation 
functions of these operators are related to the tree-level scattering amplitudes of string 
theory in highly curved AdS. The limit in which one of the R-charges is large is similar to 



1 An earlier attempt in this direction [58] related some of their building blocks to certain open 
spin chains. We depart from this standpoint and use only closed spin chains, which are perhaps 
more natural from a string theory perspective. 
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a plane wave limit. In this regime the corresponding string amplitudes for massless states 
have been computed in ||60|| . A naive extension of this calculation to massive states leads 
to an apparent puzzle: the three-point function involving one massive string state and two 
massless states vanishes identically. However, one of the examples we will discuss in §5 
leads to a nonvanishing result. Though we will not attempt it here, it would be interesting 
to clarify the relation between the correlation functions we discuss here and string theory 



computations, perhaps along the lines of pi 



2. J\f = 4 Yang-Mills and Spin Chains: the Operator Map 

The computation of planar one-loop anomalous dimensions of scalar operators in 
arbitrary representations of the R-symmetry group is conveniently encoded in the diago- 
nalization of the Hamiltonian of the Heisenberg spin chain with SO(6) symmetry 0, 

A L 

A= I^S^+ 1+2/ ^+ 1 -2^+i , (2.1) 

i=l 

where P, K and / are, respectively, the permutation, trace and identity operators acting 
on the tensor product of nearest neighbor spins i and i + 1 and L is the length of the chain. 

The spin chain description also provides a tool for the computation of various nu- 
merical coefficients appearing in correlation functions. Such techniques are particularly 
efficient in situations when the position dependence of the correlation function is known. 
The three-point functions of operators with definite conformal weights fall in this category. 
For operators of conformal dimension A/, Am and Ak the position dependence on the 
plane is fixed by (super) conformal invariance to 

/-iIMK ( \ \ 

(0 I (x)0 M (y)0 K (z))- 



(2.2) 

The OPE coefficients C IMK (X) are in general a series in the 't Hooft parameter 



C IMK (\)=c I MK (l + \c[ MK +0(\ 2 )) . (2.3) 



The one-loop correction c[ MK was partly analyzed in |58[] . Before discussing this correction 



however, it is important to have an efficient scheme for the computation of the leading 
coefficient Cq MK . 
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It turns out that the leading OPE coefficients have a simple description: since they 
are given by free field contractions, they can be identified with the matrix element of one 
of the three operators in ( |2.2| ) with the in- and out-states corresponding to the other two 
operators 

(I|6 M |J) = (M|6 7 |J) =etc . (2.4) 



IMK 




Here O m is the representation of O m on the spin chain. There is no canonical way to 
choose which operators in ( |2.2[ ) correspond to the in- and out-states. As we will see later 
it is useful to use this freedom to simplify the calculations. For example, operators with a 
small number of fields have simple expressions in terms of the basic spin chain operators 
and it is usually convenient to put them in the middle. We reach the conclusion that 
to compute Cq MJ we should compute the transition induced by some arbitrary operator 
between two eigenvectors of scale transformations, in some sector closed under RG flow. 

An alternative reasoning leads to the conclusion that the leading OPE coefficients 
can be expressed in terms of the convolution of three scalar products of spin chain states. 
Indeed, given three spin chains of arbitrary lengths in some real representation it is possible 
to define an object similar to a three-string vertex which maps a tensor product of two 
states into a third one while taking into account planar free field contractions. In the case 
of extremal three-point functions this construction can be easily extended to spin chains 
in complex representations. 

The advantage of this realization of cq is that it also provides a relatively simple way 
of computing the subleading coefficients cqC\ and more generally the coefficients of the 
logarithmic corrections to the tree level three-point functions, in terms of the expectation 
values of certain local and nonlocal operators. With a similar goal, the authors of ]58 | 



introduced a map describing the splitting of a closed spin chain into two open spin chains. 
The construction we will describe at the end of this section bypasses this step by employing 
only the spin chains describing the operators whose three-point function we are interested 
in computing. 

2.1. The SU{2) Sector 

The holomorphic operators are the easiest ones to analyze in the context of our dis- 
cussion. The reason is that the dependence of the correlation function on the in- and 
out-states becomes easy to track due to R-charge conservation. This becomes particularly 
easy for operators belonging to an SU(2) sector of the theory, 

ii...i n = Tr[ ( ^ 1 ...</> 1 ™] , (2.5) 
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where ik = 1,2 and (j) 1 are holomorphic fields. 

The dilatation operator in this sector is a restriction of ( |2.1|) and acts on these oper- 
ators as the negative of the Hamiltonian of the SU{2) XXX Heisenberg spin chain, 

L 

H = J2 H i,i+i H h2+1 = l-P ht+1 = -(a+a- +l +a-a+ +l +ala* +1 ) , (2.6) 
i=i 

with periodic boundary conditions, where of = -^(aj ±iaf) and a\ are the Pauli matrices 
acting at the site i. The eigenvectors of H can be explicitly computed for short operators, 
or for very long operators, i.e. L — > oo. In either case they are linear combinations of ( |2.5| ) 
with coefficients determined by the Bethe equations, which we will review in §3. 

Since the operators (|2.5| ) are linearly related to the eigenvectors of ( |2.6| ) which are 
special configurations of spins on a chain, it follows that to compute the OPE coefficients 
we have to find a representation of some arbitrary operator in terms of spin chain variables. 
It is not hard to see that there is a qualitative difference between operators with the same 
number of holomorphic and antiholomorphic fields and those in which their number is 
different. For the first kind of operators, ordinary Feynman diagrams imply that the result 
of the contraction of the antiholomorphic fields with fields in the in-state is an operator 
with the same number of fields as the original operator and thus it is described by a state 
in the original Hilbert space. The operators of the second kind lead to a change in the 
number of fields and thus the resulting state belongs to the space of states of a chain of a 
different length. We will treat the two cases separately. 

Operators with the same number of (f) 1 and 4>j 

It is quite easy to find the action of such an operator on the eigenvectors of ( |2.6|) (or, 
more generally, on a holomorphic operator) by analyzing the Feynman rules of the gauge 
theory. Since we are interested in the planar limit of the gauge theory, let us consider the 
operator 

O^L Tv\y* ...<:>"<:>„... , (2.7) 

which is schematically represented in figure 1 using the cyclicity of the trace. 

We see that, at the planar level, this operator takes n ordered neighboring fields with 
SU(2) indices {ji,---,j n ) an d transforms them into n ordered neighboring fields with 
indices (ii, . . . ,i n ). If ( p77|) acts on an operator which does not contain the sequence of 
indices (jx, . . . ,j n ) then the result vanishes identically. On a case by case basis the spin 
chain representation of the operator ( |2.7| ) can be written quite explicitly. 
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Jn-1 



Fig. 1: Graphic representation of the operator ( |2.7| ). 

From figure 1 it is relatively clear how to do this: due to the cyclicity of the trace 
we see that the k-th field in (|2.7|) is transformed into the 2n — k + 1-th field. Depending 
on which one of (j) , (J) 2 or their antiholomorphic versions these fields are, this transition is 
realized on the spin chain by one of the SU(2) generators. For example, using conventions 
compatible with (|2.6|) and §3, we find that 



Tr^ 1 



^ 2 



(1 + *J 



Tr[(/) 2 



E^ 1 



{21 



which are the number operators for cj) 1 and 4> 2 , respectively. A more involved example is 
the one containing all operators with, say, n fields: 



iv n [ w y + w 2 <j> 2 } n + vi +1 _^ 2 ] 
j=i i=i 

L i+n 

En 



,22 



i=l j=i 

with the coefficients u related to v and w by 



(2.9) 



(2.10) 



and the site labels defined cyclically i ~ i + L. Computing the matrix element of this 
operator between some fixed in- and out-states leads to the generating function for the 
tree level coefficients Cq MJ with fixed I and J (corresponding to the in- and out-states) 
and M of fixed size, but arbitrary SU (2) structure. 
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Operators with different number of (f) 1 and <pj 

As discussed before, these operators are qualitatively different from the operators with 
the same number of holomorphic and antiholomorphic fields, because their action on the 
spin chain changes its length. Indeed, the operator 

schematically represented in figure 2, takes n ordered neighboring fields . . . ,j n ) and 
replaces them with m ordered neighboring fields (Zi, . . . , l m ) changing therefore the number 
of fields. Because of this fact, such operators cannot have an expression similar to ( |2.9|) , 
but they also must contain a map between the spaces of states of spin chains of different 
lengths. 



Fig. 2: Graphic representation of the operator ( |2.11| ). 

In the context of our discussion it is relatively easy to construct this map based on 
the behavior of the operators (|2 .11 ). The idea is to project the in state on states which 



contain the sequence (ji, . . . ,j n ), replace this sequence with (li, . . ., l m ) and sum over all 
such states. The result is that we have to compute the expectation value of 



L-l 



fe=0 
all i s = l,2 



(2.12) 

where L is the length of the chain corresponding to the in-state, the sum is over all 
i\ = 1, 2; ... ik-i = 1, 2; ik+ n +i = 1, 2; . . .ijq = 1,2; and for k > L — n the states are 
defined using the cyclicity of the trace. Each of the two states appearing in every term 
in the sum above can be obtained by acting with SU(2) generators on the out and in 
ground states of the spin chains of length L and L + m — n, respectively, and must be unit 
normalized. Here . . . 0* n ) = |0 n ) (g> . . . <E> \(f) %n ) where l^ 1 ) = |0) and \(p 2 } = cr~\0). 
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It is worth mentioning that the representation ( |2.12|) extends to operators with the 
same number of holomorphic and antiholomorphic fields ( |2.7| ) by just picking m = n. 
Depending on the context, it may be that one of the two representations ( |2.9| ) or ( |2.12[ ) is 
more convenient. Clearly, ( |2.12| ) leads to OPE coefficients involving only scalar products 
which may be better-suited for a general analysis. 

An additional useful point is that the equation ( |2.1 2|) can be conveniently written in 
terms of the shift operator along the chain, r(0), 

T(0)\(p n (p 12 ...cf) iL ) = \(p lL (p 11 ...(/> iL - 1 ) . (2.13) 

For the purpose of computing correlation functions of gauge invariant operators (repre- 
sented by eigenvectors of the spin chain transfer matrix) this is particularly useful because 
these operators are invariant under such shifts. Thus, using 

L-l 

° l A'± <— > E r(O) fe |(^\.0^)0 i "+^.^}((0 jV .^j0 in+1 ..0 ii |r(O) & (2.14) 

k=0 
all i s =l,2 

the summation over k leads only to multiplication by the length of the chain. 

2.2. The SO (6) Sector 

The discussion in the previous subsection extends, with some modification, to oper- 
ators in the SO (6) sector, i.e. all scalar operators in the theory. The main difference is 
that, unlike in the 577(2) case, the fields transform in a real representation and therefore 
all fields in the operator whose matrix element we want to compute can have a nontrivial 
contraction with the fields in the in- and out-states. The only constraint comes from its 
trace structure. 

As before, given an in- and out-state with kj and kj fields respectively, an operator 
with kM fields will have a nontrivial matrix element between the given in- and out-states 
if 

ki = m + n kj=m + p kM = n + p . (2-15) 

Similarly to the 577(2) case, it is useful to make a distinction between the cases n = p 
(or kj = kj) and n ^ p (or kj ^ kj). In the first case we will have again two apparently 
different expressions O m which may be useful in different contexts. 
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If n = p the operator O m can be schematically represented as in figure 1, except that 
now we have to sum over all possible ways of choosing n neighboring fields to be contracted 
with the in-state. The operator 

O m = Tr^ 1 . . . </> 12 ™] (2.16) 

is represented on a spin chain by 

L 2n L 2n 

E £(^)? +1 - fe • • • (^-1)%::^ = E <°y E(^)K" +1 - fe ■ ■ • (^te-^o)-*. 
j=i fe=i ^=1 k=i 



(2.17) 



where (-E'n)i are the generators of the general linear group 



{{E n )% = 5{5] (2.18) 

with the index n labeling the site and all indices are defined cyclically 

% + L ~ % i k+2n - ik ■ ( 2 - 19 ) 

If n 7^ p the operator O m can be schematically represented as in figure 2, with the 
same provision that we have to sum over all possibilities of choosing n neighboring fields 
to be contracted with the fields of the in-state. The analog of equation ( |2.14j) is 

L-l 

n+p i y ^ T (0) k \((p ls + 1 . . 4> h +v) 4> ik + n +K . . . <j) lL \T(0) k 

(2.20) 



fc=0 
all ».=!.. .6 



where, as before, all indices are cyclically defined 

i + L ~ i h+n+p ^ l 8 ■ (2-21) 

It is certainly possible to extend this construction to operators containing fermions. 
Perhaps the simplest such sector is the one containing two scalars and one fermion (the 
gaugino in an M = 1 language), which has a manifest SU (1\2) symmetry. The extension 
to the full theory (the full SU(4\4) spin chain) is complicated by the fact that the Feynman 
diagram interpretation of the dilatation operator is not immediately obvious. We will not 
attempt to discuss this here and leave it for future work. 
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An important observation is that the technique used to construct the matrix elements 
of operators with n ^ p and different number of holomorphic and antiholomorphic fields can 
be also used to describe the mixing of single- and double-trace operators. Indeed, equations 
analogous to (|2.12| ) and (|2.20| ) represent maps between the space of states of a single chain 



and the tensor product of the spaces of states of two different chains of different lengths. 
Following the logic outlined before, the leading OPE coefficient describing the transition 
between a single-trace and a double-trace operator is given by the matrix element of 



E w a+w ■ ■ ■ <p 1s+l > ) modL2 ® w a+i ■ ■ ■ <p 1s+l ^ ) mod Ll • • ■ (<k +1 • • • k +Li ) ■ 



'IL I 



3 = 

all i n = 1 . . . 6 



(2.22) 

between an eigenstate of a chain of length L and the tensor product of two states belonging 
to chains of lengths L\ and L2 = L — L\, respectively. In the same spirit, the leading OPE 
coefficient for three operators described by spin chains in a real representation (such as 
the vector representation of SO (6)) is given by the expectation value of 



-1L 2 -1L 3 -1 



E E E E r(o)-K^...^)(^...^))® 

si=0 S2=0 S3=0 all i,j 

r(0)^|(^ . . . . .^))((0 Wl . . .0i tl )(0 im+1 • • -^ 2 )|r(0) 



(2.23) 



•S3 
£3 



between an eigenstate of a chain of length L\ and the tensor product of two states belonging 
to chains of lengths L2 and L3 obeying a relation similar to ( |2.15| ). 

2.3. The One-Loop Correction 

The form given by (|2.23| ) from the tree level OPE coefficient lends itself to higher 
order calculations as well. For the sake of simplicity we will exemplify this at the one-loop 
level, but the extension to higher loops is conceptually relatively clear. 

There are two different types of Feynman diagrams contributing to a one-loop three- 
point function: 1) the loop involves fields belonging to only two of the three operators and 
2) the four fields belong to three operators. Since only F-terms contribute nontrivially to 
three-point functions ||62|| , in both cases the R-symmetry index structure is given by the 
one-loop Hamiltonian (|2.1| ), except that the sites it acts on may belong to different chains 

H : V} a) ® Vj b) -> ® V L {d) H = H , (2.24) 
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where a, b, c, d label the chain while i, j, k, I label the sites. The different types of di- 
agrams have a different position dependence. The first type leads to the same position 
dependence as that of the one-loop corrections to two-point functions, In \xi — Xj\ 2 , while 
the second type leads to the position dependence typical for one-loop corrections to three- 
point functions, In \x{ — Xj\ 2 \xi — Xk\ 2 /\xj — Xk\ 2 ■ Identifying the labels of the positions X{ 
and the length of the operators Li, the coefficient of In |xi — X2I 2 in a generic three-point 
function is given by the expectation value of 

Li-l L 2 -l L 3 -l 

E E E E T (°y L \ ® T(o)?Hi:^\(r . . . . . . ^x^wx . . . 

si=0 S2=0 S3=0 all i,j 



|(<^ . . . <^V' fc+1 ■ • • • • • ^)>((^ m+1 ■ ■ • feJ^W • • • 0i, 2 )l r(0) 



(2.25) 

Similarly, the coefficient of In \xs — Xi| 2 |x3 — X2\ 2 /\xi — X2I 2 is given by the expectation 
value of 



L1-1 l 3 -i L3-1 



)l 2 



E E E E <oy L \®T(o) 

si=0 S2=0 S3=0 all 

• • • -^3-1^)1 T (°) 



■ / ■/ ■/ ■/ ■/ ■/ ■/ ■/ 

iL 1 jm + l im + ljL 2 irn.+l3L 2 



S3 



(2.26) 

The generalization to the other terms appearing in the one-loop corrections to three-point 
functions is now obvious. For the case of extremal three-point functions^ there are yet two 
more contributions due to the fact that the possible Feynman diagrams can have more 
general topologies (see e.g. ]65|]). These new contributions are similar to ( |2.25|) and ( |2.26| ) 
and can be easily inferred from Feynman diagrams. 

We will now proceed to describe the calculation of matrix elements of some of the spin 
chain operators constructed in this section. For the sake of simplicity we concentrate on the 
577(2) sector, but some of the formulae extend to larger collections of operators. Though 
we will mostly be concerned with the representation ( |2.9| ) of gauge theory operators on 
the spin chain, the general expressions which we will review in §4 are all the necessary 
ingredients in the calculation of the expectation values of (|2.14|) , fl2.20f) , (|2.22|) , ( |2.23|) , 
( ggp and ( fT26D as well. 



2 It would be interesting to extend this analysis to higher point functions, where extremality is 
likely to simplify the calculation. In the context of AdS/CFT extremal and near-extremal n-point 



functions of chiral operators have been discussed in |63| 1 and [p4| respectively. 
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3. Algebraic Bethe Ansatz: a Quick Review 

The SU(2) sectors of M = 4 SYM and some of its deformations are described by the 
spin 1/2 XXX and XXZ spin chains. The Hamiltonian of the XXZ chain is 



H = Yjtf°»x + + A («fi - !)) 



(3.1) 



i=i 



where erf, a = x,y, z are the Pauli matrices acting at the site i and for A = 1 we recover 
the XXX chain. This model can be solved by the Bethe Ansatz which we will now briefly 
review following |59|1, [|66 . 



The starting point is the 12-matrix which satisfies the Yang-Baxter equations. In our 
cases, the relevant solutions are 



R(X) = -i 



fa(X) \ 

/3(A) 7 (A) 

7 (A) /3(A) 

\ a(X)J 



(3.2) 



where A is the spectral parameter or rapidity, while 



a(X) = i 



for the XXZ chain and 



. sinh(A + 2irf) 



sinh 2ir] 



/3(A) =i 



sinh A 
sinh 2ir\ 



7 (A) = i 



(3.3) 



a(X) = X + i 



/3(A) = A , 7 (A) =i 



(3.4) 



for the XXX chain. The i?-matrix is normalized such that -R(O) is the permutation operator. 
The monodromy matrix is constructed out of the R matrix as 



1 ai;fci..iiM\ T>a L k L T>t>Lk 



l n h2h ' l l\)= I ^nh—nx (A) BnL'.'.'.n\ (A) 

xn £ _x • ••>n 1 W - ^ C^-^(A) D^ L ;; k n \(X) 



riL...ni 



(3.5) 



where the indices (/cl, • • ■ , ^i) and (ni, . . . , ni) correspond to the sites of the chain while 
<il and b± correspond to an auxiliary space of the same dimension as the dimension of the 
vector space on the sites of the chain. As a consequence of the Yang-Baxter equations, the 
monodromy matrix satisfies 



R(X - fj)(T[X) ® T(p)) = (T(p) ® T(X))R(X - y) 

12 



(3.6) 



where R and ® act in the auxiliary space. Examining the various entries of this matrix 
equation it is relatively easy to find that the off-diagonal entries of T act almost as cre- 
ation and annihilation operators with respect to the diagonal entries. The relevant matrix 
elements of (|3. 6|) can be written as 



(3.7) 



[B(A),S(a*)] = [C(A),C7(/*)] = , 
[5(A), C(fi)\ = g(X, n)(D(\)A(v) - D{p)A(\)) , 
D(fi)B(X) = /(A, (jl)B(X)D{h) + g(fi, A)5(/x)5(A) , 
A( M )5(A) = /(/*, \)B(\)A(p) + g(X, fi)B^)A(X) , 

where 

/(A,M) S /(A-M) = |g^ , 9 (A,,) S9 (A-,) = -I^| . (3.8) 

A consequence of the Yang-Baxter equations is that the transfer matrix, which is the 
trace of the monodromy matrix in the auxiliary space 

r(A) = TwT(A) = A(X) + D(X) , (3.9) 

commutes with itself at arbitrary values of spectral parameters. Thus, if the Hamiltonian is 
among the Taylor coefficients of t(A), the transfer matrix is the generating functional of an 
infinite number of mutually commuting conserved charges. The Hamiltonian is constructed 
out of the transfer matrix as 



H xxz = ism2r]^-\iiTxxz(X) - Lcos2r] H xxx = ^-InTxxx(X) 



d 



dX 



XXX 



d 



\=ir] dX 



- L 

A=0 

(3.10) 



and A = cos rj. 

To construct the eigenstates of the transfer matrix and the Hamiltonian one starts 
with the (pseudo) vacuum |0) which satisfies 

A(A)|0) = o(A)|0), D(A)|0) = d(A)|0), C(A)|0)=0, (3.11) 

where 

a(A) = (-ia(X)) L , d(A) = (-iP(X)) L (3.12) 

for both the XXX and the XXZ chain. With this starting point the natural ansatz for the 
eigenvectors of the transfer matrix is 

N 

l<M{A})> = n 5(A0|0), N = 0, 1, . . . , L (3.13) 
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where the arguments A^ are determined by the requirement that (|3.13 ) be indeed an eiven- 
vector of (|3l3|). Using the fact that ( |3.13| ) is symmetric in A^ it is possible to find the action 
of the operators A(fi), D(fj,) and C(n) on |^jv({A})) 

N N N N 

A(riHB(X j )\0)=A(»)HB(\ j )\0) + J2^nB(ri J] B(Xj)\0) , 

j = l j = l n=l j=l t j^ n 

N N N N 

D(ril[B(\ J )\0) = ~A(ril[B(\ J )\0) + J2^nB(ri J] BiX^O) , (3.14) 

j = l j = l n=l j=l : jjL n 

N N N 

C(n)l[B(\ j )\0) = y £Mn II BiXjm + J^MknBin) J] S(A,)|0) 

j = l n=l j = \j^zn fc>n j = l,j^k,n 

where 

N 

Mji) = a(jj) Yl /(jWjAj), An = a(A n )sr(A n , //) JJ f(X n ,Xj) , 

N 

A(fi) =d(n)Y[f(\ j ,fi), A n = d(\ n )g(fj,, A n ) ji [ /(A is A n ) , 

M n = g(fi, \ n )a(fi)d(\ n ) Yl f(Xj, A»)/(jU, Aj) + #(A n , n)a{X n )d{n) JJ /(A i5 fi)f(\ n , Xj) , 
M kn = d(X k )a(X n )g(ii, X k )g(X n , /i)f(X n , X k ) JJ f(X j: X k )f(X n: X j ) 

j^k,n 

+ d(X n )a(X k )g(fi, X n )g(X k , /j,)f(X k , X n ) JJ /(A i? X n )f(X k , Xj) . 

(3.15) 

and 

a(A) = a(A - fc) d(A) = d(A - fe) (3.16) 

with an arbitrary constant k. This shift does not affect any of the functions / and g above 
because they depend only on differences of rapidities. The standard choices, which make 
the equations most symmetric, are k = i/2 and k = irj for the XXX and XXZ chains 
respectively. 

It thus follows that |^jv({A})) is an eigenvector of ( |3.9| ) if Aj satisfy the Bethe equations 
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and the corresponding eigenvalues of the transfer matrix are 

r(M)|*iv({A})) = (A( M )+A( / .))|* J v({A})) . (3.18) 

An additional important constraint comes from the requirement that the eigenvectors 
correspond to gauge invariant operators. This translates into the requirement that they be 
invariant under shifts along the chain. Using the fact that -R(O) is the permutation operator, 
it follows that r(0) generates such shifts. Thus, the last constraint on the eigenvectors is 

A(|) + A(|) = l . (3.19) 

4. Computation of Correlation Functions 

To compute the various expectation values and scalar products derived in §2 it is 
necessary to find a realization of the spin operators in terms of the entries of the monodromy 
matrix. Such a representation makes it possible to use the relations ( |3.7[ ) and ( |3.14| ) to 
perform computations. 

4-1. Inverse Scattering Method 

The inverse scattering method expresses the local spin operators in terms of the matrix 
elements of the transfer matrix. The idea of this method is relatively easy to understand. 
It is based on the fact that there exists a value for the rapidity A for which the i?-matrix 
becomes the permutation operator. Using its properties we can see that 

T(0) = PqlPol-i ...Poi= PqxPilPil-x • ■ ■ Pi 2 = P01P12P23 ■ • • Pl-il = P01 r(0) . 

(4.1) 

Since P 01 = Yll=i a$ <8> af, we find |67| that 

a+ = C(0)r(0) L " 1 trf = S(0)r(0) L " 1 

i(l + ^)=A(0)r(0) L - 1 i(l-a 1 2 )= J D(0)r(0) L - 1 , ^ 

where we used the fact that translation by the length of the chain acts as the identity 
operator!. From here it is easy to find the spin operators at the site i, using the shift 
operator r(0): 

o^r^yo^r 1 ■ (4-3) 



3 Note that this expression holds for all models based on an i?-matrix satisfying R(0) = P 
and monodromy matrix constructed as in (p.5|). Thus, the invariant meaning of the arguments of 



A, B, C, D and r in equation ( |4.2| ) is the value of the spectral parameter for which the i?-matrix 
becomes the permutation operator. It will occasionally be convenient to shift the rapidity, as in 
equation ( 3.16j) . 
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It is however not hard to see that, using these expressions, the equations ( |3.14|) lead to 
apparently singular results if we act on some state with two or more operators of the same 
rapidity. Indeed, both /(0) and gr(0) are singular, and going back to equation ( |3.6| ) does not 
yield the necessary commutation relations. Using the explicit form of the spin operators 
one can easily convince oneself that the scalar products are completely finite. Thus, for 
computational purposes, it is necessary to regularize (|4.2| ). A convenient regularization is 
provided by the anisotropic chains, which are based on the same .R-matrix as the isotropic 
one, except that the i?-matrices building the monodromy matrix are evaluated at site- 
dependent shifted rapidities 

T(X) = R ol (X-Cl)Rol-i(X-Cl-i)...Roi(X-(i) , (4.4) 

and the vacuum state ( |3.11| ) and the vacuum energies ( |3.12|) are modified appropriately. 
Though the arguments are more involved ] 68|| , the expressions for the spin operators 



in terms of the entries of the monodromy matrix are quite similar to Ql2| ) and 

i-1 L 



a=l a=i+l 
i-1 L 

*t = {]l(A + D)(Za)}C(Z i ){ n (A + D)(C a )} 

a=l a=i+l 

l{l + *) = {f[{A + D){Z a )}A(Z i ){ f[ (A + D)(Z a )} 

a=l a=i+l 

kl-<r?) = ill(A + D)(Za)}D(Z i ){ f[ (A + D)(£ a )} 



(4.5) 



a=l a=i+l 



According to [67] similar formulae hold for superalgebras, which generalizes our discussion 



to some sectors of gauge theory containing fermionic operators. 
4-2. Scalar Product 

As we have seen in §2, the calculation of correlation functions of gauge invariant 
operators reduces at the level of the spin chain to the calculation of the scalar products of 
states constructed out of B and C operators 

N N 

^ = <0| C(Af ) n i?(Af ) |0> , (4.6) 

J=l k=l 
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where Xj are kept arbitrary (i.e. not necessarily satisfying the Bethe equations). There is 
a large amount of literature devoted to the computation of such inner products, part of it 
described in detail in [59]. One of the main results is that ( ^.6| ) is given by 



Sn = II S(A?, X°)g(X B , A?) sign(P c )sign(P B ) J] h(Xf B , Af c ) \{ h(Xf c , A 



H MAf C , Af c ) n h(\f B , Af B ) det(M^) det(M^) 



Lfe 



(4.7) 

where P c and P B are the permutations {Af c , . . . , A^ c , Af c , . . . , X BC i n } of {Af , . . . , A^} 
and {Af B , . . . , Af B , Af B , . . . , Xjf_ n } of {Af , . . . , A*} respectively, and 

(M$E) jk = t{Xf B , X^)d{X^)a{Xf B ), t(A, p) = 9 4k^, h(ji, A) = (U 



A) 



This result is derived using the commutation relations (|3.7|) as well as various properties 
of rational functions. When the in-state is an energy eigenstate (i.e. {X B } is a solution of 
the Bethe equations ([3.171) ) (|4.6| ) simplifies considerably and takes the following form |38| 

S N{{x c h {x b }) _ IL d(Xf) Ua d(X B ) det H({Xf }, {Af }) 



n i>fc v(Af-A?)ni< m ^-Af; 



(4.9) 



where 



<p(v) 



¥>(A£ - Af ) 



o(Af) 



d(A c 



; n ^(Af - Af +77o) - n ^ - x * - »») 



(4.10) 



The function 99(A) equals A and sinhA for the XXX and XXZ chains, respectively, while 
r/o equals i and 2ir\ for the XXX and XXZ chains, respectively. 

The entries of the matrix H have potential singularities if some rapidity Af of the 
out-state approaches one of the rapidities A^ of the in-state. This apparent problem is 
however cured by the fact that, in this limit, the square bracket in ( |4.10| ) turns out to 
vanish due to the Bethe equations (|3.17| ). This type of scalar product is sufficient for the 
calculation of three-point function coefficients and will be useful in the examples we will 
discuss in 55. In the limit in which the two states become identical, (]A. 91) becomes the 



square of the norm of the Bethe eigenstates given by the Gaudin formula proven in 

(4.11) 



N N 

n c(Xj) n B(A fc ) io> = n d^d^ n /(a oj a 6 ) det $'({a» 
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where $' is an iV x N matrix with the elements 



d 



dX, 



log 



N 

n 



/(A„, AO 



<W t ijU/(A*,A»; 



(4.12) 



The explicit form of these matrix elements is not too complicated and simplifies even more 
in the thermodynamic limit (N, L — > oo with ^ fixed) as well as in the long chain limit 
(L —>■ oo with iV fixed). The explicit form of (|4.12|) is 



ab 



Sab 



d 



(A„)-~(A a ) + ^^(A a -A, 

k=£a 



+ K(X a -X t 



where 



(4.13) 



(4.14) 



It turns out that K is the kernel of the integral equation for the density of rapidities in the 
thermodynamic limit. After using the Bethe equations in the same limit, the coefficient of 
the Kronecker 5 becomes the density of rapidities. Combining everything, <£' is 



Kb(W) 



<W(A a ) + ^K(X a - A b ) 



(4.15) 



The ratio N/L enters here through the normalization of p. 

In the long chain limit the result is even simpler because the solutions of the Bethe 
equations are of the order of the length of the chain, A a = L/2-Kn a . Thus all terms 



depending on K in (|4.13| ) have subleading contributions and it is not necessary to solve 
the equation for p. The matrix $' becomes diagonal, with entries given by the first two 
terms in ( [4.13| ) and its determinant is therefore 



det$'({A}) = n 



a=l 



a ' / > N d' ,y, s 

A a ) + - A a 

a a 



(4.16) 



Each term in the product is of order of the inverse length of the chain. At first sight, 
this dependence as well as the extra factors in ( [4.11| ) seem to be in contradiction with 
the known results concerning long operators. Let us therefore briefly discuss the relation 
between Bethe eigenstates and gauge theory operators. 

Since we picked a normalized vacuum state |0), the corresponding gauge theory oper- 
ator must also be normalized. We therefore have 

1 



10} 



Vl 



Tr[(0 



l\Li 



(4.17) 
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If creation operators are present, the solution to the apparent puzzle implied by ( |4.16| ) is 
provided by the cyclicity of the trace and the fact that the action of the creation operators 
B(X) brings in additional rapidity-dependent factors. The relation between (|3.13| ) and the 
non-normalized gauge theory operators is, roughly, 



N 



i^M{A})>=n s ( A *)i°> 



i=l 



N 



Y[d(\, 



i=l 



' N 

n 



O 



non — norm 



VI 



(4.18) 



where (9 non - norm j s? as usual, a linear combination of traces of products of (L — N) (ft 1 
and iV (j) 2 fields while yL in the denominator comes from the normalization of the vacuum 
state ( |4.17| ). Thus, Y\*i=i d(Xi) cancels the similar factor in the norm Q4.11|) of |\Pjv({A})), 
while the other factor contributes L~ N . Thus, combining (|4.18 ) and (|4.16| ) we find that 
the normalized BMN operators are related to the non-normalized ones as 



1 



O 



non — norm 



(4.19) 



Using the cyclicity of the trace we can always place a <p 2 field as the first field inside the 
trace; this reduces the number of terms in (9 non - norm by a factor of L and simply multiplies 
the resulting operator by L. We therefore find 



l*iv({A})> 



o, 



non — norm 



L (N-l)/2 



non— norm 



^Tr[^ p (^,^ 2 ,{A})] 



(4.20) 

which is in agreement with the known gauge theory results. Here K p denotes some mono- 
mial in the fields with a A-dependent coefficient. It is certainly possible to do a more 
detailed analysis and keep also the subleading terms. We will see this explicitly in an 
example in the next section. 



5. Examples 

Let us illustrate the method of computing correlation functions with some examples. 
All examples will involve states created by two creation operators, so we first normalize 
them. For a chain of length L the corresponding solution of the Bethe equations ( |3.17| ) 
can be explicitly found for the XXX chain 

Ai = -A 2 = -ictg-^y , (5.1) 
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where A^ are the shifted rapidities ( |3.16| ). Then equation ( |4.11 ) implies that the normalized 
state is 

\^ oim (n)) =Af 2 B(X 1 )B(-X 1 )\0) 

Ai /. 2 1V~ L A L ~ l Tin . 2fL _ n Tin ( 5 - 2 ) 

M 2 = ; A? + - = ; COS Sin 2(L l) . 

y/L{L-l) V V ^L(L-l) L-l L-l 

Extracting the operators B out of the monodromy matrix and acting with them on the 
vacuum state leads to an explicit expression for the (normalized) gauge theory operator 
associated to the Bethe state |\E , 2° rm ( n )) 



L — 2 

O n = — L= £ cos ™j> + 1} Tr^ 2 ^)^ 2 ^)^ 2 -] , (5.3) 



s=0 



in agreement with known results [ 70 ] . For other spin chains some simple normalized gauge 
theory operators corresponding to Bethe states have been constructed in |14| . Calculations 
similar to the one leading to ( |5.3| ) yield (after solving the Bethe equations) all finite length 
operators in the SU(2) sector with definite one-loop anomalous dimensions. 

Example One 

A simple example of the type in which two operators have the same length is provided 
by a three-point function of a generator of the Cartan subalgebra of the i?-symmetry group 

Tr^i - cp 2 fo] (5.4) 

and two other operators. According to the discussion in §2, the tree-level coefficient of this 
correlation function is given by the following expectation value: 

L 

(n oim (M)iE^i^ orm ({A})) = ^(n orm (M)i^in orm ({A})) . (5.5) 

Using the solution ( f4.2|) of the inverse scattering problem this matrix element reduces to the 
matrix element of A(0) — -D(O) between the in- and out-states above. In this computation 
it is not necessary to use the regularization provided by the inhomogeneous chain, because 
there are no potential singularities from contractions. This matrix element was computed 
in |68[ and the result is the one expected on physical grounds 

L 

<n orm (M)i E^i^ orm (W)> = ( L - 2N Wi x } - M) » ( 5 - 6 ) 



i=l 
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where <5({A} — {//}) vanishes unless the set {A} is the same as the set {fx}. Indeed, the 
operator ( |5.4j) counts the difference between the number of (j) 1 and (f> 2 fields. Since the 
action of B replaces (p 1 with </> 2 , the result should be the difference between the number of 
(p 1 fields in the vacuum and twice the number of creation operators. Furthermore, since the 



Hamiltonian commutes with (|5.4[) , the latter does not mix the energy eigenstates, which 
explains the delta function. 
Example Two 

A more complicated example is the extremal three-point function 

C&[0f- a ](a;)Tr[g](y)O 2 W> , (5.7) 



where O2 is the normalized operator corresponding to \^2° Tm ( n )}- Using equation ( |2.14| ) 
it is not hard to find that the tree level normalization of ( |5.7| ) is 

2 

c =L L -A0\<p ll ...<p lL - 2 )(^2K...^ L . 2 m° Tm (n)) L , (5.8) 

il...i_L_2 = l 

where the intermediate states are considered to be normalizedB, the indices on the vacuum 
states represent the length of the chain describing them and the overall factor of L comes 
from the sum over the powers of the shift operator in ( |2.14j) . The scalar product between 
the out vacuum state and the intermediate state forces the intermediate state to also be 
the vacuum state \0) L _ 2 . Thus, only one term in the sum survives 

c o = LM' 2 (fafa$ 1 ...$ 1 \B(\ 1 )B(-\ 1 )\0) Ii , (5.9) 



where A/2 is the normalization constant of the in-state ( |5.2| ). Using the solution ( |4.2| ) of 
the inverse scattering problem and the fact that ( |5.1| ) represents shifted rapidities, the 
remaining scalar product reduces to! 

c = L Af 2 (0| Ca )CU ) B(X 1 )B(-X 1 )\0) . (5.10) 



4 In fact, it is not hard to see from §2 that the intermediate states are already unit normalized. 

5 As discussed before, the arguments of the C operators above must be the values of the 
rapidities for which the i?-matrix becomes the permutation operator. For the shifted rapidities 
this is Aq = §. 
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Using (|3.14| ) and the regularization provided by the anisotropic chain or the general ex- 
pression for the scalar product between an energy eigenstate and an arbitrary state (|4.9|), 
we find 

L-l 



(0|C(i)C(i)S(Ai)B(-Ai)|0> = 2 A 2 + 



1 



2 3-2L sin 2(l-L) 



Tin 
L-l 



(5.11) 



Combining this with the normalization of the in-state (|5.2|) leads to the three-point function 
coefficient 

/ L 4Ai / L nn 

Similar coefficients have been computed for some operators dual to string modes in 
plane wave background. The equation Q5.12 ) is more general since it is valid for operators 
of arbitrary length. 

As a side remark, it is not hard to see using the analysis of §2 that the same tree-level 
three-point function coefficient arises, up to some numerical coefficient, if we replace Trj^ 2 ,] 
with Tr[(0 1 ) fc (^|] for any k < L — 4 and also appropriately adjust the length of the chain 
describing the out-state. 




two-point function type 



three-point function type 





three-point function type; topology special to extremal correlators 



Fig. 3: Various diagrams contributing to fl5.7f) at one- loop level 



Using the analog of the equations Q2.25 ) and (|2.26|) for the SU(2) sector it is not hard 
to find the coefficients of the one- loop corrections to (|5.7| ). The relevant gauge theory 
Feynman diagrams, grouped as discussed in §2, are listed in figure 3. The white dots 
correspond to fa while black dots correspond to fa. Some of the contributions are partic- 
ularly simple. Since both Tr^^ -2 ] and Tr^ 2 ,] are BPS operators, they are annihilated by 
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the Hamiltonian of the spin chain and the contribution of the analog of equation ( |2.25| ) 
(first two diagrams in figure 3) vanishes. The contribution of the analog of ( |2.26|) does not 
vanish and it is a sum of two terms. The simplest contribution (third diagram in figure 3) 
comes from the identity operator in the Hamiltonian and contributes cq to the coefficient 
of In \z — x\ 2 \z — y\ 2 /\x — y\ 2 . A similar contribution comes from the interactions typical 
to extremal three-point functions (fifth and sixth diagrams above). 

The final contribution (fourth diagram in figure 3) comes from the permutation oper- 
ator and is proportional to 



C£> = (o|C(|)(A(|) +J D(|))C(f)|^r m H) 



(5.13) 



This matrix element is a special case of a spin-spin correlation function and measures the 
probability of having two down spins with an up spin between them. The computation of 
this scalar product is relatively easy. As before, we need to regularize it by shifting the 
arguments of the operators in ( |5.13| ) at arbitrary positions. It is not hard to see that the 
nontrivial normalization factor of the regularized state is irrelevant, so we will not compute 
it explicitly. Definingi 



s(o = (o\c(oc(mr im (n)) 



A lv /TT4A|(±-z£) 



L-1 



it 



L-1 



(5.14) 



y/L{L-\){\\ - e) 

which can be easily computed using (|4.9|) , the scalar product ( |5.13| ) is given by 



4 1 ) = (a(!) + rf(t))S(t) + 



(5.15) 



S 2 



Using the explicit form of a, d and S(£) we find that the last one-loop three-point function 
coefficient is given by 



t>P — Co 



1-4 sin^ 



Tin 



1 



L 3Tcn 
cos - 



1 



1 



(5.16) 



Collecting the contributions of all diagrams in figure 3 yields the one-loop correction to 
the three-point function (|5.7|). 



Clearly, S(- 



c . 
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